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ABSTRACT 


A study was conducted of the amplitude modulation of an 
acoustic signal forward Scattered from a wind-generated, 
model ocean surface. The demodulated spectrum of the surface 
scattered sound. was compared with the model ocean spectrum 
while varying the acoustic frequency and angle of incidence. 
Several methods of enhancing signal detection through a 
knowledge of the ocean spectrum were studied. The most 
promising of these shows potential gains of the order of 
three decibels when the received signal is weak and 


fluctuating greatly, as at a Lloyd Mirror minimum. 
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I. INTRODUCTION 


It has long been known that received signals from 
distant underwater sound sources vary in time. The effect 
has been described as "builds and fades" of the signal. 
Though this effect is automatically interpreted by the 
human ear, which listens when the sound is present and tunes 
out when it disappears in the noise, little or no advantage 
of the phenomenon has been taken by the extensive signal 
processing equipment avallable today. Though nominally far 
superior to the ear in response-time and frequency resolution, 
this equipment often fails to achieve even parity with a 
human observer when processing a time varying, or modulsted, 
Signal. This is due primarily to the electronic equipment's 
narrow bandwidth and the fact that it integrates continuously 
rather than selectively. 

It was the objective of е research to investigate some 
of the parameters affecting this phenomenon of a modulated 
acoustic signal with particular attention to determining, for 
Surface affected signals: 


I DheMbestelabor2story modelling techniques. 


о те гетасјог од (пе шоапГаштоп frequency to the 
surface wave frequency. 

Be ine Variation of the modulation amplitude апа spectrum 
with frequency and source/receiver directivity. 


ши ре optimization of Signal processing of an amplitude 
modulated signal by digital techniques. 


Aa: 





In short an attempt was made to discover a "temporal 
filtering" scheme which could be used in concert with 
existing high speed frequency filtering techniques, such as 
the Fast Fourier transform, to enhance signal processing. 


The concept has been dubbed "EASE", Environmentally Adjusted 


Signal Enhancement. 
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IE THEORY 


The subject of scattering of sound from the ocean 
surface has been treated by several authors utilizing nearly 
as many methods. The following is a summary of Refs.1 and 2 
which point up those features of the scattered spectrum 
which are germane to understanding the relationships between 
sound fluctuations in the upper ocean and the ocean wave 
f spectrum. 

As Parkins summarizes the phenomenon: "When a mono- 
chromatic acoustic wave is incident on the surface of the 
ocean the envelope of reradiation varies in time or fades." 

Assuming that the ocean is generated by a stationary 
homogeneous, Gaussian process, and that it is smooth enough 


so that multiple scattering can be neglected, the incident 
TOt 


acoustic wave of harmonic time variation e is given by: 
> > 
= er 
"E - am 1 „Tut 
1 
with 
КУ Ка 
= бон [2x + my + nz) 


pt issan emplitude term, ® m, and n are the direction 
cosines and À is the wave length of the incident acoustic 


wave. Other factors are defined in Figure 1. 
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Z = t(x,y,t) 


INCIDENT 
ACOUSTIC 
ENERGY 





FIGURE 1. Coordinate System for Theoretical Development 


Following the development of Beckmann and Spizzichino 
[Ref. 2] which stipulates a Source free, isovelocity ocean, 
and using the Helmholtz Integral development of Stratton 


[Ref. 3], the farfield reradiated pressure may be expressed 


аз: 


14 





1(wt -KR,) R. 





К г > 
le iL 1K‘? 
sem I [Г 4- dxdy 
r 218. K, L 
with: 
EEC тук Б у. x 
е y? С ЗІ 


and other factors again as defined in Figure 1. 

The irradiated area has been assumed to be an L xL 
aperture which is several wavelengths across. 

In this research it was the spectral density of the 
reradiated pressure that was of interest. This may be 


obtained by first forming the autocovariance function, 
R(t) = «6, (x,y,t) 9 F(x,y,ttt)» 


The autocovariance function having been found, the Wiener- 
Khinchin Theorem states: 


~ 


Scu uM Mni aU ст 
= СО 
where 2 is the frequency of the reradiated (or scattered) 
sound. 
Continuing with Parkins, and assuming that the surface 
irregularities are much smaller than the aperture, the 


СТЕБЛО OT Бе spectral density of the scattered sound 


Ша 


mE 





is as given below. It should be noted at this point that 
the requirement for small surface irregularities corresponds 


to a roughness parameter, 
E = Ко = сто [соз 0- + соз 0. ] 
Z À 1 2 


as defined by Beckmann and Spizzichino, much less than one. 
o is the ocean RMS wave height, A the acoustic wave length 
and 01 and 9, the angles of incidence and reflection, 
respectively (see Section IV.A.4 for a further discussion 
or yes). 

The expression for the spectral density of the scattered 


sound is: 


-K 292 16 sin“ 


50) = Be = [5(9-ш) 5 5 
К, Ку 


2 
K 2 sin Kym 


Dac e | 2 2 2.1 l 
+ 2% K, o {[A (KK) | S[N-w-(K, tK, ) (06): 4 
+ [A'(-K,,-K,) J^5[0-u*(K, ^«k ^) *(0) 311 
Хати X y 
yip cU uemaccdleraulsoncobemsua vu 
A' the Neumann-Pierson directional wave spectrum 


B the amplitude term 


СК "ге Dirac delta 


16 





Parkins summarizes this on: "o. (this equation) 
shows the signal reradiated by a wind generated sea to have 
a specularly reflected component occurring at the incident 
frequency w and two scattered components which are Doppler 
shifted by equal amounts from о." 

The spectrum of the reradiated sound is thus seen to be 
frequency modulated by the Doppler shift term and amplitude 
modulated by the wave spectrum. It is the amplitude 
modulation by the ocean waves that is considered in this 


research. 


E 





III. EXPERIMENTAL SETUP 


This research was carried out using the Naval Postgraduate 
School OCEAN ACOUSTIC WAVE FACILITY (OAWF) with data 
collection and analysis accomplished with a hybrid computer 
setup designed by the Special Projects Group of NADC and 
affectionately referred to as OPHELEA (CCEAN PHYSICS 
ENVIRONMENT EFFECTS ANALYZER). All other equipment utilized 
was standard commercial electronic gear and the author feels 
that a detailed explanation of its functions and capabilities 


is unnecessary (see Equipment List - Section III.D). 


А. ОАМЕ 

This facility eonsists of the unique combination of a 
wind-water tunnel and an anechoic tank (see Figures 2 and 3). 
Up to five 3/4 HP centrifugal fans (three in this case) 
comprise a Sea Exciter (Figure 4) which forces air down the 
wind-water tunnel, which is about 17 meters in length, 
generating waves which empty into an anechoic tank. The air 
space in the tunnel was standardized at 14.25 cm and the 
Sroscm coe) Ome ol sene Bentiire Sbpunnel measures about 1.2 meters 
square. Observation of the waves at various points along 
the tunnel is permitted by glass observation windows shown 
in Figure 5. Waves entering the 3x3x2.9 meter tank portion 
of the facility are traet Irom Со со ОН Oya 


plexiglass restraining barrier and are terminated by a 
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"beach" to prevent reflections. The beach spans the width 
of the tank experimental area and consists of a Screen 
covered wedge filled with stainless steel millings which 
cause incident waves to dissipate rather than reflect. 

The tank itself is lined with redwood 4x4's on all sides 

and the bottom which ensures excellent sound absorption at 
the boundaries as pointed out in Ref. 4. These redwood 
beams are oriented with a corner facing toward the tank 
interiors except on the wall bordering the waves which 15 

| flat, to further increase absorption. А movable race permits 
two-dimensional positioning of transducers to within one 
centimeter. This facility models the real ocean on a Scale 
of about 50-1 and is shown in operation in Figures 6 through 


8. 


В. OPHELEA 

The OPHELEA system consists of essentially three 
components which are interfaced to provide high speed analog 
to digital data collection and digital data processing with 
somewhat slower display. The design was developed by the 
Special Projects Section of the Naval Air Development Center 
ШП соп ШШ шшш tchsPaenkertonetemputer Consultantsftliner, 
O Warmimstverssrennsyıvaniasss Thesthree сотпропепоз б аге вһомп 
in Figure 9 and described below. 

1. Interdata Model 70 Computer 

This minicomputer is modularly constructed and 


features LSI and other state-of-the-art computer technology. 
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It is a 16 Bit half-word oriented digital computer that is 
FORTRAN programmable and has sufficient memory (32 thousand 
8-bit bytes) capacity to process in excess of 8000 data 
points with an average size program. In this research it 
was programmed to perform a Fast Fourier Transform on input 
data and output spectral density vs. frequency. Some 
statistical analysis was also performed. This minicomputer 
occupies the upper half of the large cabinet shown in 
Figure 9. 

2. Phoenix Analog to Digital Converters Model ADC 712 

Two ADC 712 converters may be used separately or 

simultaneously. Each A/D converter is a high speed device 
capable of encoding up to +10 volt input signals into 
12 bits (11 bits + sign) of data, providing resolution to 
one part in 4095 (accuracy of 5mv in 20v). The maximum rate 
of conversion is 6.5 microseconds and it has performed 
reliably at a sampling rate of up to 300kHz. The A/D is 
driven by a "command to convert" signal consisting of a 
5 volt positive pulse of at least 200 microseconds duration 
which is supplied by an external oscillator. Measurement 
is accomplished by successive comparisons against an internal 
reference voltage after which the value is converted to a 
digital number. The A/D converters are located in part of 


the lower half of the OPHELEA cabinet. 
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3. Texas Instruments Silent 700 Electronic 
Data Terminal Model 733 


The TI 733 consists of a keyboard which is used as 
a programming and input/output control device, a printer 
er outputs data or answers queries of the computer on 
heat sensitive paper and a play back and record section 
which is used to read in computer programs and record/ 
playback data. In this research, due to storage limitations, 
data were sometimes taken A/D, read out to tape cassettes, 
then reread in for processing after the appropriate program 
had been loaded. A capability of interfacing with the 
IBM-360/67 CP/CMS system, allowing direct access to the NPS 
computer facility, was also utilized in processing some 


statistical wave height information. 


С. SEQUENCE OF OPERATION 
The sequence of operation used most often in this 
research was as follows: 
| 1. An analog signal (varying sound pressure, instan- 
taneous wave height, etc.) was input to the A/D converter 
which, being triggered at twice the Nyauist rate (a value 
found to result in minimal allasing effects), converted it 
Od CICA dala 
2. Those data were then either: 
a. directly processed by the appropriate program 
Unle овез сае ке computer, Ог, 
ШО read out to tape саззебтез Tor later analysis 


as described above. 
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3. Reduced data and desired results were displayed on 
the paper readout of the terminal. Figure 10 is a pictorial 


description of this process. 


D. EQUIPMENT LIST 

In order to avoid needless repetition of the lengthy 
manufacturers! names for electronic equipment throughout 
this work, a number of abbreviations will be used instead. 
These short titles will be employed for the most part in 
block diagrams and schematics but may also appear in the 
text. The abbreviations, with the corresponding full 


equipment title/description are as listed below: 


ABBREVIATION MANUFACTURERS TITLE 

CORRELATION COMPUTER Princeton Applied Research 
Model 101 Correlation Function 
Computer 

DEE POVERSSWEERTZ Hewlett-Packard Model 721A DC 
Power Supply 

G.R. OSCILLATOR General Radio Model 1312 Decade 
Oscillator 

HP PRECISION AMP Hewlett-Packard Model 463A 
Precision Amplifier 

HP POWER AMP Hewlett-Packard Model 467A 
Power Amplifier 

KH FILTER-3342 Khron-Hite Model 3342 Filter 

KH FILTER-3550 Khron-Hite Model 3350 Filter 

®' SCOPE Tektronix, Inc. Model 545B 
Oscilloscope 
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ABBREVIATION 


PAR PREAMP 


WAVETEK 142 
WAVETEK 144 
X-Y RECORDER 


MANUFACTURERS TITLE 

Priiszeton Applied Research Model 113 
PreAmplifier (Note: HFRO and LFRO 
refer to the High/Low frequency roll 
off feature of the preamp — а 
filtering option) 

Wavetek Model 142 HF VCG Generator 
Wavetek Model 144 HF VCG Generator 


Varian Aerograph Model F-100 Strip 
Chart Recorder 
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IV. EXPERIMENTAL PROCEDURES AND RESULTS 


A. OAWF CHARACTERISTICS 

| In order both to ensure the validity of the research to 
follow and to gain an insight into the significant parameters 
affecting the scattering of sound from a real or model ocean, 
several preliminary experiments were conducted. In the 
sections that follow each will be fully discussed with setup 
diagrams, procedural explanations and results included. 

1. Frequency Spectrum 
In determining the frequency Spectrum of the ocean 

wave heights, a highly linear probe designed at the Civil 
Engineering Department of Stanford University and pictured 
in Figures ll and 12 was utilized. The probe is capacitive 
in nature, with a covered wire as the sensor. The wire 
serves as one electrode, the water as the other and the 
insulation is the dielectric of the capacitive element. 
The probe is excited by a 2.4 kHz sinusoidal carrier 
(amplitude 4.5 volts peak to peak). When the probe is 
immersed in the water wave system the output across 4 
capacitive bridge (Figure 13) is amplitude modulated by 
the’ pasising waves! rlsing and falling around the wire and 
thereby changing the effective length of the capacitor. 
This AM signal is demodulated by the detector shown in 
Figure 14 and the resultant time varying DC signal, a replica 
of the instantaneous wave height at one point, is frequency 


analyzed by the OPHELEA system. 
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Figure 11. Stanford 
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Figure 1H. 
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As noted above frequency analysis was accomplished 
by the OPHELEA system, using a program called PING IC, a 
copy of which is included in Appendix A. In this case the 
program called for taking 256 A/D samples at a sampling rate 
of 20 Hz and computing a 128-point spectral density versus 
frequency output. These parameters provide a frequency 
resolution of greater than .08 Hz. 

When taking A/D samples at sampling frequencies in 
excess of 20 kHz, a UE oscillator which is half-wave 
rectified by computer circuitry was used to provide the 
"command-to-convert" pulse train. At low frequencies however 
the sinusoidal pulse duration is too great and erroneous 
multiple triggering of the A/D Conver ver eres] Comme Laermer Oe 
the WAVETEK 144, with a positive pulse waveform and with its 
asymmetry control set for minimum pulse duration, was employed 
for the low sampling frequencies required for finding spectra 
of wave and demodulated scattered sound. A block diagram 
of the experimental setup is shown in Figure 15 and the 
normalized spectral density plot is shown in Figure 16. 

Ths plot iS thesayerage Ofeforty runs conducted as outilined 
above (i.e., 40-128 point FFT analyses). 

As can be seen from the figure, the predominant ОАМЕ 
ocean frequency is about 2.5Hz as compared with a typical 
Acmucalmececanmmreaqueney On tne order of 0.182. Though neither 
ocean is monochromatic, the OAWF spectrum is fairly narrow. 


Its bandwidth is about 0.6Hz at the half power points and 
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the spectral density tends to fall off less rapidly toward 
higher frequencies than does the actual ocean Spectrum 
[Ref. 5]. This is due to the greater contribution of 
capillary waves in the OAWF waves. The irregularity at 
about 2.4Hz was felt at first to be merely a coincidental 
artifact which would disappear with additional averaging. 
However, though it smoothed some, it has persisted. The 
source of the irregularity may have been lack of adequate 
fetch for wave growth to an equilibrium spectrum. It is 
noted that a similar effect was observed in the work of 
Fowler and Scheibel [Ref. 6] at short fetch. 
горар тиеу Density Function 

As Kinsman has noted [Ref. 7] the probability density 
function associated with wind driven ocean wave heights is 
nearly Gaussian. In order to verify that the OAWF sea met 
this criterion, wave height samples were collected using the 
setup shown in Figure 15 of the previous section and analysis 
was accomplished using the standard IBM-360/67 Library 
program "HISTF". This process was accomplished as follows: 

a. Eight groups of 102! wave height samples were 
taken from the previously described wave probe using the 
OPHELEA system and a program modestly entitled PERK 14 and 
were stored in the minicomputer memory. 

b. Тһе data were then read out to tape cassettes in 
a revised format compatible with the IBM-360 CP/CMS system. 

c. Data were converted to punch cards using the 


SENS” system and analyzed using "HISTF". 
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The "HISTF" program outputs various statistical 
parameters and plots a histogram of the input data. Each of 
the eight data groups was analyzed and the results averaged. 
A graph of the resulting average histogram, as well as a 
comparable normalized Gaussian PDF, is shown in Figure 17. 

As can be Seen, the histogram distribution is slightly skewed 
and displays negative kurtosis. Though not exactly agreeing 
with any one of Kinsman's seas, the variations are slight as 


is shown by the following comparison: 


Kinsman - record 072 OAWF - 3 Fans 
Skewness +, 092 –. 0245 


Kurtosis = 0: -, 3549 


The mean is also Seen to be non-zero but is less 

than 0.1 centimeter. 
Э. Autocorrelation Function 

Ihe autocorrelation function of the ocean wave 
heights is characteristically that of band limited noise and 
takes the shape of an attenuated cosine [Ref. 8]. In finding 
pnexuPnalatrow runctionwotbthe OAWF three fan ocean, 
essentially the same setup was used as in the previous two 
sections. However the output from the wave probe, 
eharaceterizime che surfaces was terminated, not in the 
OPHELEA system as before, but in a PAR CORRELATOR. This 
device is capable of computing either auto- or cross- 


ВоВе Ome runections*and displaying the result on an 
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oscilloscope or a recorder. In this case an X-Y RECORDER 
was used. A diagram of the experimental setup is shown in 
Figure 18. 
An average surface wave height correlation function 
о пеш та ВУ 0o UO SECCONOS was arrived әс ру 
superimposing the eleven runs printed on the X-Y RECORDER 
and shown in Figure 19. The resultant average is shown in 
Figure 20 and can be seen to follow the expected form. The 
rate of "attenuation" is initially somewhat less than 
measured by Kinsman, with the OAWF correlation at the first 
peak down to 0.7 compared to 0.50 for Kinsman's ocean data. 
The second peak correlation conversely has dropped to about 
0.3 while Kinsman's C(t) remains at about 0.4 after time 
Hagst of two ocean periods. 
The correlation funetion loses its consistency after 

two ocean periods as may be seen from Figures 19 and 21. 
Figure 21 shows the function computed with a precomputatlon 
period = T = 2.0 seconds set into the PAR CORRELATOR. With 
рініс setting the correlation function is computed for 
2.0 to 4.0 seconds or in effect the function of Figure 19 
ДИ Convinueu for another two seconds. For this longer time 
Баса 15 сітізеліз Еос ееси пу repeatability in Successive 
correlation functions. 

While the OAWF values are somewhat different from those 
ef the Kinsman ocean; the similarity in shape is obvious. 


This validation of the ocean model will be of importance 
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Figure 18. Experimental Setup for Determining the OAWI 
Surface Wave Autocorrelation Function 
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later in this work when comparisons are made with the 


modulation envelope of scattered sound. 


4. The Roughness Parameter 


The previously mentioned roughness parameter, є , 


is defined by Beckmann and Spizzichino [Ref. 2] as: 


Ve = ате [cos 04 + cos 9 ] 
with o = ocean RMS wave height 


à = acoustic wave length 


9, = angles of incidence and reflection, respectively 
(measured from the normal to the ocean surface). 


При! сисази ту seen that a "smooth" condition (Ve << 1) сап 


be achieved through use of low acoustic frequencies or large 


angles of incidence (shallow grazing angles) as well as with 


the obvious small RMS surface wave heights. 


The converse is 


of course also true for "rough" conditions (Jg >> JT). 


Ме was the factor that was used to model the OAWF ocean to 


the actual ocean in later research and an example of typical 


вова тон № Боби отајптпо умејатлпе thefsame 


МЕ = .076 is given below: 


Wind 

RMS wave height 
Acoustic Frequency 
Source/Revr Depth 


Source/Revr Separation 
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OAWF 


3 Fans 
ШЕ СП 
КТЕ 
10 см 
150 cm 


OCEAN 


DS 
54 cm 
SOO THES 
200 m 
7000 м 





As noted above Mg. Hs the) same for both sets of 
parameters and thus the typical ocean conditions noted in 
the right column can be reproduced if the values in the left 
column are employed in the OAWF. 

B. COMPARISON OF WAVE HEIGHT AND DEMODULATED SCATTERED 

SOUND SPECTRA AND CORRELATION FUNCTIONS 

As noted in the Theory section the frequency spectrum of 
the amplitude modulation of the scattered sound is directly 
related to the frequency spectrum of the wind driven waves. 
Hence it is important to compare these frequency spectra. 
Further it is of interest to compare autocorrelation functions 
in order to gain an insight into signal processing possibili- 
ties. In addition to these two comparisons the cross 
correlation function was computed. 

1. Frequency Spectra Comparison 

Though some trials described later in this work were 
conducted with a directional source, the principal study 
used an omnidirectional source and receiver, which were felt 
to best model the submarine/sonobuoy scenario. The source 
WaS & locally-constructed transducer made Lon a 2 ке са. 
Glennite Ceramic sphere which was covered with several coats 
of neopreme. This source transducer vill be referred to 
hereafter as "OMNI", The receiver was a standard Atlantic 
Research LC-10 hydrophone and its response was immediately 
Emp gc ОЕ ООШ ОТ Ога бй FET preamp powered Dy 


ШЕШ POWER SUPPLY Both the LC-10/NUS combination 
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and the OMNI, as well as a directional piston type transducer 
(F-27) used later, are shown in Figure 22. 

As the spectra of the OAWF sea waves had already 
been computed from an average of forty runs, only the 
average sound amplitude modulation spectrum remained to be 
determined. To accomplish this the scattered AM sound was 
demodulated by the detector shown in Figure 14 and the 
resulting time-varying DC signal sent to the OPHELEA systen. 
A photograph of the wave forms of the AM scattered Sound, 
the demodulated scattered sound and the instantaneous wave 
height are shown in Figures 23 and 24. A complete schematic 
of the experimental setup is shown in Figure 25. 

Forty runs were conducted with 256 A/D samples of 
the demodulated acoustic 20 kHz Signal being taken for each 
run at a sampling frequency of 20 Hz, again using program 
PING 1C. As before, a 128-point spectral density versus 
frequency readout was calculated for each run and the 
resulting forty spectra averaged and normalized for comparisons 
with the previously determined average wave spectrum. These 
two averaged, normalized spectra are plotted in Figure 26 
which shows the obvious similarity. The demodulated sound 
spectrum has a shape that is remarkably similar to its 
Counterpart eneuch its peaks ame about 0.2 Hz less then the 
predominant wave frequencies. This difference is presumably 
ШІС our surface -draive ve locitys| Ref. 9]. The 3 dB bandwidth 
of the demodulated sound spectrum is about 0.5 Hz as compared 


ею Неро» Бле wave spectrum. 


48 








Hydrophone (LC-10) and Transducers Used 


in this Research 


Flgure 22, 





ВА ве № 


Amplitude Modulated 
опа op) and 
Instantaneous Wave 
пете и 





Figure 24. 


Demodulated AM Sound (top) 
and Instantaneous Wave 
Heicht 








OMNI LC-10 


RE | HP | | GR-1312 
r ا‎ POWER AMP < = OSCILLATOR 


AMP=1K ( | 
= Hz PAR 

KH FILTER | 

! ха 


$ ОО пай 





| NUS PREAMP AND 
wai DC POWER 
| SUPPLY (12 v) 










DETECTOR | 
(See Fig. M) — 





vas 


“ЗАМР-.02К 


KH FILTER 


3 5 5 0 ! Шен ты 


LF 


I LFR0=DC 





PAR | X-Y 


CORRELATOR Lemna RECORDER 





Figure 25. Experimental Setup for Determining the Spectrum 
of Demodulated Scattered sound 
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These. frequency spectra, though exhibiting some minor 
differences, are considered’*tG*2ssxssentially the same and ` 
confirm the theoretically expected result. 

20, Autocorrelation Function Comparison 

The same procedure was used in determining the 
Gemodulated scattered sound autocorrelation function as in 
the case of the wave height version. The output of the 
previous setup was terminated in the PAR CORRELATOR instead 
of OPHELEA and recording was done on the X-Y RECORDER as 
before. Again eleven runs were conducted for delays 
UE to 2.0 seconds. These recordings are shown in 
Figure 27. Ап average autocorrelation function was then 
fitted to this record. A plot of both records, wave height 
and demodulated scattered sound, is shown in Figure 28. 
Though the ordinate at T = O0 is slightly different for the 
two functions due to idiosyncracies of the PAR CORRELATOR 
and X-Y RECORDER, it is apparent that the results are nearly 
identical for the two cases. Further, it is evident that 
s Y они оостгед ас тони ро смила ав евнеха съв вот three 
"peaks", that is, over two ocean periods. The demodulated 
sound correlation function also becomes confused for delays 
Брсагет than 2.0 seconds’, as) is"the case with the wave height, 
as seen in Figure 29. It should be noted that the vertical 
Scalewnskigure 29 is@expanded Over that of Figure 21 and по 


infenence of greater eörrelation should be drawn. 
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3. Cross Correlation Function 
Using the two Setups .descrihed in the autocorrelation 

section above the two wave forms were cross correlated, again 
uSing the PAR CORRELATOR and X-Y RECORDER. The same sort of 
averaging process was then performed on the record shown in 
Figure 30 and ylelded an average cross correlation function 
which is shown in Figure 31. The t = 0 ordinate varied 
considerably between positive and negative values of 
correlation and the choice of a near zero initial value was 
purely arbitrary. Again three peaks are seen to dominate, 
at about t = 0.8, 1.25 and 1.7 seconds. This is to be 
expected from the cross correlation function of two signals 
which both exhibit significant correlation over two ocean 
cycles (or have three "peaks" in their autocorrelation 
functions). 
C. EFFECT OF FREQUENCY AND ANGLE OF INCIDENCE 

ON SCATTERED SOUND MODULATION 

To determine whether the modulation frequency of the 
scattered sound remained essentially constant with changing 
acoustic frequency and varylng angle of incidence, a series 
of trials was conducted with different combinations of these 
two parameters. It was felt at this stage of the research 
that there might be some advantage in using a directional 


hydrophone to make the detection. 
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пилане пе Frequeney and Angle-of itneidence 


In these tests a directional transducer, a USRD 
Type F-27, was mounted in a trainable bracket which permitted 
angles of incidence, measured from the surface normal, of 
from 0° to 90°. The F-27 is shown in Figure 22, The 
scattered sound was received by an omnidirectional LC-10 
hydrophone which was positioned at the same depth as the 
F-27 and 2.75 meters away. To achieve different angles of 
incidence while retaining specular scatter the F-27 and 
LC-10 were positioned at depths of 0.99, 0.49 and 0.27 meters 
which equate to angles of incidence of approximately 58°, 
70° and 79° respectively. Prior to each trial the F-27 was 
elevated and depressed, with the surface quiescent, until 
the maximum response ensured that the observation was of 
specular scatter. Ten runs, using program PING 1C (Appendix 
A), were conducted at each position with 256 A/D samples 
being taken at a sampling rate of 20 Hz. A 128-point spectral 
density versus frequency output was generated as in the 
determination of previous spectra. Again the WAVETEK 
oscillator provided the pulse train for proper A/D converter 
trigecering at low*frequencies. "А diagram of this Setup is 
shown in Figure 32 and a table of run parameters and 
particular equipment settings is found in Table II. 

After the above 80 runs were completed, the ten 
spectra from each anfle/freaueney combination were averaged 


and normalized to permit comparison. The results were then 
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plotted as indicated in Table I below with results compared 
for constant angle of incidence and variable frequency and 


vice versa. 


TABLE I 


SUMMARY OF VARIABLE ANGLE OF 
INCIDENCE AND FREQUENCY TRIALS 


Figure Acoustic Frequency Angle of Incidence 
33 30,40 kHz 79° 
34 20,30,40 kHz 70° 
35 20,30,40 kHz 598 
36 20 kHz 7O ee 
37 30 kHz Том 5,59 
38 HO kHz | A DOS 


As may be seen from Figures 33 through 38 the 
frequency spectral peaks differ by less than 0.3 Hz in all 
cases except Figure 35. The somewhat wider spread, about 
0.5 Hz, between peaks is thought to be due to greater 
contributions trom the bade lobes at that angle or incidence, 
Even this) difference, thoughelarcer Mis comparable to the 
3 dB bandwidth of the wave frequency spectrum. 

Based on this work, 1t was concluded that the 
modulation spectra of the scattered sound received by a 
Docneanydeophome saremsimilar for variations in acoustic 
| ir l o Roer a о or incidence, It is felt that the 
differences are due to inadequate averaging and/or the 


dependence on side lobes. 
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2. Modulation Strength with Changing Receiver Angle 


The object of this set of, trials was to determine if 
the modulation effect could be increased by training a 
ЯН кессіспеі recelver to a preferred angle of reflection with 
a point source. Although acoustic reciprocity would be 
expected to hold, in this case control would be at the option 
of the listener. Specifically it was felt that possibly if 
the receiver hydrophone were directed normal to the modulating 
ocean Surface, the relative modulation of the received 
scattered acoustic signal would be increased. The experi- 
mental setup was essentially зе before; however, now the 
F-27 was the receiver and the previously referred to OMNI was 
Еге SOuUrce WIth bot аса јаерсп Of 1.0 meters. sine’ three 
fan OAWF sea was again the environment with the water level 
standardized as before. А setup diagram is shown in Figure 
39. Ten runs were conducted at each of seven angles of 
incidence ranging from direct path, through specular to 
normal incidence. The acoustic frequency was 30 kHz and 
amplifier and filter settings were kept constant for all 
runs so that valid comparisons could be made. Table III 
Enowswssummarysom theseastrials. 

Баси сЕ cher seventevarunszinzsTtaple Lil consisted of 
taking 256 A/D samples at a 20 Hz sampling frequency (using 
program PINC 1C) as before; then averaging the spectral 
denservsinlormation for each angular setting. The plot of 
these averare spectra is presented in Figure 40 which shows 
lopespcerral density plotted vemsus scattered sound 


modulation frequency. 
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TABLE ЕТ 


SUMMARY OF VARIABLE ANGLE. RECEIVER TRIALS 


Run Numbers Angle of Incidence 
61-70 90° (Direct) 

71-80 58° (Specular) 
81-90 150 

91-100 202 

101-110 0° «(Normal 
111-120 255 

121-130 152 


Note: Again the angle of incidence is measured 
normaleroschezsurfsceh 


All levels have been plotted without adjustment 
So that the amplitude of the modulation at different angles 
of incidence can be seen Simultaneously. Summarizing: 

a. The greatest modulation strength was observed 
at 6 = 58°, that is the specular direction. 

b. Modulation amplitude declined regularly as the 
angle of receiver is changed in either direction from the 
SPECULAT 

с. Although the incoherent component (modulated 
amplitude) may possibly be increased relative to the 
coherent component, the modulation amplitude is down 30 dB 
(re specular) with the receiver trained normal to the sea 
surface (0 = 0%), 

Knowing that the amplitude modulation is greatest at 


or near the specular direction might provide a means of 
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Figure 39, 
of Chanringe Receiver Angle on Modulation Strength 
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arriving at a rough fix or target probability area. Assuming 
that target bearing could be provided by a directional hydro- 
phone and that the local sound velocity profile is known, the 
approximate range could be found by applying ray path 
techniques once the angle of greatest modulation (specular 
direction) is found. This Scheme would of course require 

a hydrophone system which is directional in both azimuth 


and elevation. 
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V. ENHANCING THE TIME-VARYING SIGNAL 


At this stage, with reasonable assurance that the OAWF 
sea accurately modeled the actual ocean and that the modula- 
tion spectrum of the scattered sound was essentially that of 
the surface waves, efforts were begun toward improving signal 
processing using the model and the information gained to date. 
If the ocean wave spectrum were of a single fixed frequency 
the sound would fluctuate at that same single fixed frequency 
and sound signal enhancement would then be easy. It would be 
necessary merely to identify a loud point in the received 
signal and then to listen at later intervals equal to the 
ocean wave period. The ocean is regrettably not this simple 
as is obvious from the respective frequency spectra and 
therefore alternative or modified processing means had to be 


identified and implemented. 


A. PROGRAMS PERK 1A AND PERK 1B 

Mo berin tLomexplores@the problem, a pair of computer 
programs modestly named for this.rookie acoustician were 
Geveloped which would accomplish the following: 

1. PERK 1A orders the computer to direct the A/D 
converter to take a group of samples (say 64), hereafter 
сап ея а “BLOCK” ада then tom pause”. 

2. The "pause" is triggered by the computer taking its 
last ordered sample and PERK 1A then shifts the computer 
Man “Do Boop" which has am index proportional to 


the desired time delay. 





3. When the "Do Loop" has completed its preset number 
of operations (and the delay or pause has been completed) the 
program signals the A/D to take another BLOCK of samples. 

This is followed by another "pause", another BLOCK, etc. 

The net result is a series of BLOCKS of A/D samples 
Caken at specified Intervals alone =the modulated signal 
which, when a total of 8192 data points have been taken (the 
present storage limit of OPHELEA), are read out to a tape 
cassette. This data tape is then put aside as the second 
program, further modestly entitled PERK IB, is loaded, 

PERK 1B reads the time series data of the previously 
recorded sample BLOCKS back into memory destructively and 
performs an FMI spectral analysis on each BLOCK. The program 
may then be directed to output all of’ the complete spectral 
analyses, or just the major spectral peak in each analysis. 
Thus a sequence of successive sound spectral density peaks 
is avallable, each separated by a fixed time delay, the 
"pause", The duration of the "pause" is selected to provide 
the desired number of spectral densities of predetermined 
frequency resolution during an "ocean period". An "ocean 
period" is here defined as the reciprocal of the peak 
frequency in the ocean spectrum. The effect of the modulation 
caused by the ocean is vividly illustrated by Figure 41, which 
depicts the peak sound spectral density versus time; and which 
shows a swing of about 18 dB between maxima and minima of 
Sienal strength. Coptes of PERK IA and PERK 1B may be found 


in Appendix A. 
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In this work, the BLOCKS almost invariably contained 
64 A/D samples and the "pause" was selected to provide eight 
BLOCKS per ocean period. The ‘following is a sample delay 
time calculation for a 20 kHz acoustic signal, modulated by 
a 2,5 Hz ocean (ocean period 0.4 sec.) and sampled at 80 kHz. 

For 64 Samples/BLOCK at 80 х 103 samples/sec.: 
sampling time/BLOCK = 0.8 msec. 

For 8 BLOCKS/ocean period: sampling time/ocean period = 
6.4 msec. Therefore, for one ocean period, total delay time 
= 400 msec - 6.4 msec = 393.6 msec. 

Delay between BLOCKS = 393.6/8 = 49.2 msec. 


Total data collected in 16 ocean periods 


(93 samples y (8 BLOCKS у (26 осеап periods, 
BLOCK ocean period kun 


samples 
8192 run 


For the purposes of this research eight BLOCKS per ocean 
period were considered adequate coverage. Each of the eight 
BLOCKS was said to be at a different PHASE lettered A through 
H. Every eighth BLOCK would be at the same phase. Again if 
the ocean Were monochromatic and the spectral density peaks 
in each PHASE averaged, one average, corresponding to the 
strongest part Of the signal, would predominate and tne other 
А ош Lassereampditude Would bewedistributed sinusoidally 
around Lu. 

PERK IA and PERK 1B are highly versatile programs with a 
number of options available as Summarized in the following 


table (Table IV). 
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TABBER V 


SUMMARY OF PERK 1A AND PERK 1B PARAMETER OPTIONS 


Total number of A/D data points Up 1.096192 

Number of data points/BLOCK from 16 to 512 in powers 
of 2 (2N) 

Delay between BLOCKS 0.01 to 999.9 msec with 


accuracy to 0.01 msec. 


Number of BLOCKS/ocean period 2 to 16 in powers of 2 
= number of PHASES 


ї 


The program also provides a unique number for each run, 
based on day, month and run sequence of that day (i.e., 
#15101 is the first run conducted on 15 October), and also 
prints out input parameters such as sampling frequency, prior 
to displaying the output. 

It should be noted at this point that PERK 1A/1B 
were tested on a known sinusoid of 2.5 Hz and found 
satisfactory. In that test, groups of 8 A/D samples separated 
by 0.4 seconds were taken at 80 kHz and plotted versus time. 
No visible discrepancy existed over 2! sinusoidal periods. 

torre capitulate sethesprogramns! PEER LA “and PERK IB were 
used to take BLOCKS of 64 A/D samples, properly spaced 
temporally to provide eight PHASES (A-H) per ocean period, 
ver sixteen ocean periods, This equates to 128 BLOCKS or 
16 BLOCKS per PHASE. Samples were taken at four times the 


acoustic frequency, or twice the Nyquist rate, to diminish 


| 





the effects of aliasing. “A pictorial т Ton ce 
this signal processing scheme is shown in Figure 42. 
B. EFFECT OF SIGNAL STRENGTH ON MODULATION — 

LLOYD MIRROR MODELING 

One of the most obvious times when signal processing 
enhancement is important is when the received acoustic 
Signal is weak and may be fluctuating in and out of the 
background noise. In order to provide a common basis for 
comparing varying signal strengths, it was decided to employ 
different aspects of the Lloyd mirror phenomenon to observe 
varying strengths of the output signal. 

1. Brief Review of Lloyd Mirror Effect 

As every school-boy acoustician knows, the Lloyd 

Mirror Effect is caused by constructive and destructive 
interference between the direct and smooth surface reflected 
paths taken by an acoustic Signal. For a mirror surface within 
Critical Range ва - Eo the effect is one of pressure 
peaks and valleys, with the peaks falling off as 1/r and the 
valleys going to zero (h = receiver depth, d = source depth, 
АХ = acoustic wavelength). Beyond the critical range the 
pressure declines as l/r*. The OAWF tank is capable of 
measurements both within and beyond the critical range for 
frequencies in the vicinity of 13 kHz. Another means of 
Visualizing this effec: is busreplacing the reflected rays 
by an image source at I as also shown in Figure 43. 

In the case of a rough surface the peaks caused by 


relnforcement are not as great and the coherent scattering 
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Figure 43. Lloyd Mirror Effect 
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term, given by the first Dirac delta term in the expression 
for the scattered sound spectrum (Theory section), predominates. 
Further the valleys caused by cancellation are now minima 
instead of zero’ as the incoherent or modulated term in the 
S(Q) expression becomes important. This effect is shown in 
Figure 44 which contrasts the rough and smooth situation. 

The setup used in the following three trials was essen- 
tlally the same as in all signal processing runs. The OMNI 
was the source and the LC-10 served as receiver (both at a 
depth of 10 cm.), with a 3 Fan OAWF sea, set to the standard 
water level, modulating the transmitted 20 kHz acoustic 
signal. A/D samples were taken with program PERK 1A at a 
sampling rate of 80.0 kHz and analysis was performed by 
PERK 1B. A setup diagram is shown in Figure 45. The rack 
on which the LC-10 was mounted was moved back and forth with 
the water surface quiescent (fans off) until a Lloyd mirror 
maximum and adjacent minimum were identified. A third 
position was then chosen arbitrarily between the two extremes 
and designated "Intermediate". Runs were conducted at each 
of these three locations and a Summary is presented in Table V 
below. 

From the output of PERK 1B it was possible to compare the 
signal fluctuations at the three positions. At the Lloyd 
mirror maximum, corresponding to a strong signal (larger 
constant amplitude coherent and smaller variable incoherent 


contribution from the surface scattered sound), the maximum 





TABLE V 


TRIALS CONDUCTED AT VARPO S~- le Ese MIRROR POSITIONS 


Run No. Lloyd Samples/ Delay No Sas, Total Ocean 
Mirror BLOCK PHASES Periods Covered 
10101 Min 64 7 G9 8 16 
msec 
10102 Max 64 (assumed 8 16 
ocean 
Гуеа. 
10.03 Inter- 64 2,578 Hz) 8 16 
mediate 


spectral density peak exceeded the overall average by about 
2.1 dB and the minimum peak fell about 2.0 below the average. 
The total range of amplitudes was about four decibels. 
Applying the same comparison technigue to the other two runs 
yields a swing of about 8.5 dB at the "Intermediate" position 
and almost 19 dB at the position of the weak signal where 
there is a maximum fluctuation of the amplitude. 
тара вас 0115 model fluctuations qdo exist 
regardless of signal strength and that they are of greatest 
importance when the signal is weak. This further points up 
the importance of adapting knowledge of the sound modulation 
to enhancing a signal processing method as a means of 
acquiring and/or retaining low strength signals in a noisy 


or hishly variable environment. 
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Figure 45. Experimental Setup for Determining the Effect 
on Spectral Density Variation with Time at 
Three Lloyd Mirror Positions 













ШЕН "ВЕ ПШИЕНАЗЕ АВТОР 

As previously mentioned. if the. scattered sound were 
modulated at a single fixed frequency there would be a PHASE 
which, when once correctly chosen, would always yield the 
strongest average signal. The modulation frequency is of 
course not monochromatic as shown earlier in Figure 41 and 
in fact its spectrum has about a 0.5 Hz half-power bandwidth 
(Figure mo. 

This not withstanding it was felt both instructive and 
possibly rewarding to compare the relative strengths of the 
eight PHASES. Thus an annex was appended to PERK 1B which 
averaged all of the spectral density peaks in each PHASE 
(every eighth peak in this case) and printed out a PHASE 
average as well as the overall average of all peaks. These 
values for 8 runs are shown in Figures 46(a~h), which depict 
PHASE average spectral density plotted versus PHASE with the 
dashed line indicating the overall average. All of the runs 
are seen to, in fact, follow a rough sinusoidal pattern 
particularly Figures 46(c) and 46(g) but others such as 
Figure 46(a) appear to be at double the ocean frequency. 

The Figures 46(a-h) show spectral density on a linear scale. 

If the "BEST PHASE", or highest PHASE average is compared 
to the overall average, the gain is in the range of 0.3 dB 
(for bhael Lloyd mirrorfmax=strones signal case) to 0.5 dB 
(for the weak E which is hardly anything to write 


home about. 
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In summary the "BEST PHASE" method applied innocently 
over 16 ocean periods is seen to offer some small advantage. 
However, because of the broad band modulation of the scattered 
sound modulation and the drift of the best phase over many 


Ocean periods the gain realized is minimal. 


D. "MONDAY MORNING SIGNAL PROCESSING" METHOD 

The next approach began with hand plotting the spectral 
density versus time information presented by PERK 1B and 
examining it to see if any combinations would become apparent. 
Samples of the plots are included in Figures 41, 47 and 48 
and show the fluctuations of spectral density with time for 
the Lloyd mirror minimum, "intermediate" and maximum, 
respectively. 

To analyze this displayed information a clear plastic 
mask was constructed with a sinusoid of the same frequency 
as the assumed peak OAWF ocean frequency printed on it. The 
mask was then moved along the time axis of the spectral 
density versus time plot to determine how many peaks or near 
peaks corresponded to the ideal sinusoid peaks. Typically 
ron аб least MICE бо fives would Line up prior 
to an irregularity or "glitch" which would require shifting 
phase to a subsequent group. This was not unexpected as the 
auto- and cross-correlation functions discussed earlier in 
Бо ти па: сасва коса corie lation for at least Two or 


three ocean periods. 
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The above procedure was felt to simulate a signal 
processor which would identify a strong signal and then 
sample later, only at intervals which are multiples of the 
ocean period, three to five times to gain signal enhancement 
ОСИ СОО ЕЦБ И Ие сла оп The еси отисне signal 
processing in hindsight — hence the Monday Morning label — 
are interesting and are illustrated below: 

Referring to Figure 41, the Lloyd mirror weak signal 
(large fluctuations) case; if the BLOCKS indicated by arrows 


are chosen and averaged as follows: 


2 


Group BLOCKS M SP Ave. 
l DIN CNN 44,70 
e SIS OS e ppc 
3 Dior nO ore ТИМ Шә? 12.55 


These groups of BLOCKS correspond to identifying initial 
strong signal (i.e. BLOCKS 5,39 or 98), then sampling at the 
ocean frequency for only three to five peaks further. The 
payofi becomes epparent when these Monday Morning averares 
are compared to the overall Ae of 19.6 which yields 
palm о 1 апарзьз, "ваветою groups 1, 2, and 3, 
respectively. 

Gains in the case of the modeling strong signal (weak 
fluctuations) ranged from +1.0 to +1. dB and those for the 
Lloyd Mirror "Intermediate" response were in the +0.4 to 


91.1 ав range. The interested reader is invited to verify 


9! 











these values, Or try to improve on them with other 
combinations, by marking a ruler or like object in intervals 
of 0.4 seconds and sliding it along Figures 41, 47 and 48. 
The advantage of Monday Morning Signal Processing, is 
Obvious in the weak signal environment (where it is most 
needed, of course) and is at least worthwhile, though not 
ie ие о-в оп Ретро 1 спав Шлевацсротї will 
consider modifications to this technique and suggestions for 


real time analysis in the final section. 


E. SUCCESSIVE AVERAGING 
In an attempt to validate statistically the "seat of 
the pants" notion that three to five ocean periods is the 
maximum extent to which sampling at multiples of the ocean 
period can be profitably employed, a system of averaging was 
developed and appended to PERK 1B (which became PERK IC). 
The spectral density peaks were first divided into eight 
PHASES (A through H) of sixteen peaks each. These PHASE 
groups were then blindly averaged in ensemble group as follows: 


fei ahem ASE В for example, 








Average by two's zz , а A B16+B1 
T est +В3+ Eq 
Average by three's M I Ве = B4 n B16 z B2 


etc. 


Up to Average by eights 





These averages were then compared with the overall 
average and the number that exceeded it noted. Rather than 
confirming the initial hypothesis, however, ‘the results were 
contradictory; indicating that groups of three to five were 
most profitable in some PHASES while averages of five to 
eight were most beneficial in others. In fact one PHASE 
indicated that it made little difference how many were 
' averaged — that all combinations were equally profitable. 

The author remains convinced however, that, though it 
has not been statistically proven, when compared with blind 
averaging, the Monday Morning technique demonstrates that 
averaging more than three to five peaks in a group is 
counterproductive. 

In summary the blind Best Phase method provides oniy a 
minimal gain of about half a decibel when applied to a weak 
signal while the Mondary Morning technique provides nearly 


а +3 ав enhancement of the signal processing for very weak 


Sean sa mp lerne ges Uc dm Gost nice #708 Tl ve7öoceanTzeycieck 
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VI. FUTURE RESEARCH, CRYSTAL BALL GAZING AND OTHER MUSINGS 


ie author feels That Lie research described" in tHe 
foregoing pages has only seratched the surface of a viable 
means of acoustic signal enhancement. Some basic parameters 
ЕЛІНЕ сатре Шо Ene modulation вон Расе stattered Sound 
have been identified or confirmed and in some cases loosely 
quantized. These ideas may be of immediate use to signal 
processors, however additional research is felt to be 
indicated in the following two areas: 

A. Further investigation of the phenomenon for additional 
clues which might increase its value/utility and prove more 
rigorously some notions presented herein. 

The initial concept of identifying a signal peak (or loud 
point), then sampling later at multiples of the ocean period, 
which at this point might be charitably described as 
simplistic, could be more effective in a modified form. If 
a peak were identified and the next sampling interval were 
ОКО И the anterval T2¿AT” where T is the Ocean period 
and AT is related in some unknown way to the bandwidth of 
the ocean spectrum and/or the rate of "attenuation" of its 
autocorrelation function, greater improvement might ke 
realized. Of course the interval of sampling would increase 
В ооо: рез ог АТ if the starting point were not 
reset at a peak after cach group of samples was taken and 
the peak identified. (See Figure 49.) It is felt that the 


ИО ПЕТ о ес ТПІтІР Of three to five sampling intervals 
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Figure 49, 
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inproving on the 


Interval not reset to signal peak; 
next sampling interval must be 
+2AT to cover all possibilities. 
Then +3AT, etc. until continuous 
integration results 


Interval reset to signal peak 
and next peak will now be 
within Т+АТ 


"Best Phase" Processing Method 





before the diminishing returns point is reached may also 
apply here. 

It might be even more profitable to take samples at a 
rate of several BLOCKS per ocean period over some observation 
interval (greater than the ocean period) as done in this 
research and then to compute, sort and rank the spectral 
density peaks by magnitude. Once so arranged the following 
options would be available: 


l. use only the highest peak in each observation interval 


2. use only the highest peak in each ocean period and 
average over the observation interval 


3. -use the five (or any other arbitrary number) highest 
peaks observed and average as in 2.). | 


1. use only those peaks which exceed the overall average, 
which would also have to be computed, by 3 dB ( or 
employ some lesser standard, such as 2 or 3, if none 
duds 

As an aside it should be noted that when speaking of the 
"highest spectral peak" it is implicit that the highest for 
each frequency observed is computed, compared and utilized. 
This research was conducted with a single acoustic frequency 
but future trials might include several tranducers operating 
at various frequencies and signal output levels. A BLOCK of 
more than 64 samples may be required in this case to provide 
the desired frequency resolution. 

Another possibility might be to continuously process for 
mMOStmen the loud, half ef thefocean period (5 seconds or 
so in the actual ocean), then process and use the received 
frequency spectral information to image correlate during 


the quieter half cycle. Again the three to five period 
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ООО А EMIL Ss feft собресгайт, ThisSsimerchodlmight 
provide even better results st eruthssccobtained by active 
sonar Image correlating as the changing doppler effects and 
a Cendant correlacion mismatches would net De a factor over 
such short time spans. 

B. The second major area of suggested study would be to 
spectrally ae plot, and employ the Monday Morning Signal 
Processing Technique on actual passively-recorded submarine 
tapes to see if gains comparable to this research are realized. 
Additionally the modulation frequeney Spectrum and autocorre- 
lation funetions could be determined and compared with both 
the locally observed wind/sea data and the results of this 
work. 

вла ни спете не ne Со бю с СО he. fact Спар Sound 
scattered from the ocean surface iS amplitude modulated at 
the basic surface wave frequencies. The ability to estimate 
the ocean wave frequency and spectrum is of course vital to 
taking advantage of this phenomenon and has not been addressed 
pO TNS DOLCE The@author feels that ballpark estimates 
(ош решаг- туеа atl on station through û knowledge of fetch 
and surface wind speed and duration and assuming a Pierson- 
Moskowitz wave height spectrum. Further it is felt that 
meteorological organizations Such as Fleet Numerical Weather 
LEN от дос ел роса тв r inrormatcion if tasked со do so, 
БОЗИ Иол Processors would be provided with up-to-date, 


tactically-useful data. 
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оон mechanisms such as internal waves and tem- 
perature microstructure also eontribute to fluctuations, the 
surface scattering is of particular importance for near- 
surface sources and/or receivers. What remains in the author's 
opinion is to find the means to most effectively use the 
various ocean modulation effects to enhance signal processing. 
It is felt that clever signal processors, if properly in- 
formed, could devise numerous techniques far more subtle wand 
comprehensive than the rather simplistic methods suggested 
in this and the foregoing sections and that even greater 
gains might be realized. It should be emphasized that use 
of these techniques is not contended to be the miracle cure 
for the Anti-Submarine Warfare problem. The author does feel 
however that ocean modulation phenomena, if used to best 
advantage, will provide part of the answer required to cope 


ШЫЛ Ош сек next generation targets. 
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